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Abstract 



A new measurement of the inclusive production cross section for pp —> tt is per- 
formed at a center-of-mass energy of 7 TeV using data collected by the CMS exper- 
iment at the LHC. The analysis uses a data sample corresponding to an integrated 
luminosity of 36 pb^^, and is based on the final state with one isolated, high trans- 
verse momentum muon or electron, missing transverse energy, and hadronic jets. 
The tt content of the selected events is enhanced by requiring the presence of at 
least one jet consistent with b-quark hadronization. The measured cross section is 
150 ± 9 (stat.) ± 17 (syst.) ± 6 (lumi.) pb and is in agreement with higher-order QCD 
calculations. The combination of this measurement with a previous CMS result based 
on dileptons gives 154 ± 17 (stat. + syst.) ± 6 (lumi.) pb. 
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1 Introduction 

The top quark was first observed in proton-antiproton collisions at ^/s = 1.8 TeV at the Fermilab 
Tevatron collider ||TJ|2l. Since then its properties have been studied by the Tevatron experiments 
and found to be in agreement with the expectations of the standard model (SM) 13]. At the 
Large Hadron Collider (LHC) [4J top quark production can be studied in pp collisions at a/s = 
7 TeV, allowing extended measurements of the top quark properties. A precise measurement of 
these properties is important as top-quark production may be a background for new physics. 

At the LHC top quarks can be produced singly or in pairs. This paper focuses on the study of 
the tt final state, for which the production cross section has been calculated in next-to-leading- 
order (NLO) and approximate next-to-next-to-leading-order (NNLO) quantum chromodynam- 
ics (QCD) ISHZI- In the SM the top quark decays almost 100% of the time via the weak process 
t — > Wb. We focus on the tt decays in which one of the two W bosons decays hadronically and 
the other decays leptonically (the semileptonic channel), giving a final state containing an elec- 
tron or muon, a neutrino, and four jets, two of which come from the hadronization of b quarks. 
Here taus are detected only through their semileptonic decays to electrons and muons. The 
results are based on the analysis of a data sample corresponding to an integrated luminosity of 
36 pb"^ |8|, which was recorded by the Compact Muon Solenoid (CMS) experiment between 
April and November 2010. 

We present the results of tt cross section measurements in the muon and electron + jets chan- 
nels using the subset of the data in which at least one jet has been b tagged using a displaced 
secondary vertex algorithm. A profile likelihood method is used to fit the mass distribution of 
the identified b-decay vertices as a function of the jet and b-tag multiplicities in the event. The 
analysis is sensitive to the differences between the signal and background processes and also 
allows the simultaneous fit of the light (u, d, s) and heavy quark (b, c) contributions. The main 
systematic uncertainties are taken into account when maximizing the profile likelihood. This 
allows for the correct treatment of their correlations and the evaluation of the combined uncer- 
tainty. These results complement the CMS dilepton and kinematics-based lepton+jets analyses 
which are described elsewhere [9, 10 1. The ATLAS experiment has also measured the tt cross 
section in the dilepton and lepton+jets channels at 7 TeV IITTjI . Several cross-check analyses are 
also performed, which use different analysis techniques (a fit without the use of the profile 
likelihood, and two measurements based on simple cuts), and different b-tag algorithms (one 
based on a track impact parameter, and another that uses soft muons). 

The central feature of the CMS detector is a superconducting solenoid, 13 m in length and 6 m 
in diameter, which provides an axial magnetic field of 3.8 T. The inside of the solenoid is outfit- 
ted with various particle detection systems. Charged particle trajectories are measured by the 
silicon pixel and strip tracker, covering < (p < In in azimuth and \yj\ < 2.5, where the pseu- 
dorapidity rj is defined as rj = — ln[tan(0/2)], and 6 is the polar angle of the trajectory with 
respect to the anti-clockwise beam direction. A crystal calorimeter and a brass/ scintillator 
calorimeter surround the tracking volume and provide high resolution energy and direction 
measurements of electrons, photons, and hadronic jets. Muons are measured in gas-ionization 
detectors embedded in the steel return yoke outside the solenoid. The detector is nearly her- 
metic, allowing for energy balance measurements in the plane transverse to the beam direction. 
A two-tier trigger system selects the most interesting pp collision events for use in physics anal- 
ysis. A more detailed description of the CMS detector can be found elsewhere [|12|. 

We describe the data and event selection in Section 2 of this paper, followed by a brief descrip- 
tion of the modeling of the signal and background processes in Section 3. Section 4 describes 
the method used to extract the cross section from the selected events, as well as the calculation 
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of the statistical and systematic uncertainties on the result. The cross-check analyses are dis- 
cussed briefly in Section 5. Lastly in Sections 6 and 7 we present a summary of all of the CMS 
measurements and compare the results with the predictions from QCD. 

2 Event Selection 

The trigger used to select the data samples for analysis is based on the presence of at least one 
charged lepton, with either an electron or a muon with a transverse momentum, pj > p!^". Be- 
cause of increasing maximum instantaneous luminosity the minimum transverse momentum 
was varied between 9 and 15 GeV for muons and between 10 and 22 GeV for electrons in order 
to maintain a reasonable trigger rate. The same data are used both for signal selection and for 
the study of the non-top QCD multijet and W(Z)+jets backgrounds. These triggers have been 
shown to have efficiencies of 92.2 ± 0.2% and 98.2 ± 0.1% for the muon and electron channels, 
respectively, based on independent studies of Z boson decays into ^^}^^ and e^e^ pairs. No 
azimuthal nor polar angle dependence is observed. 

Muons are reconstructed using the information from the muon chambers and the silicon tracker 
and required to be consistent with the reconstructed primary vertex [13J. A kinematic se- 
lection requiring pj > 20 GeV and \f]\ < 2.1 is then used to select muon tracks for further 
analysis. Electrons are reconstructed using a combination of shower shape information and 
track/ electromagnetic-cluster matching ||T4| . A veto is applied to reject the electrons coming 
from photon conversions. To be retained for further analysis, electron candidates are required 
to have a transverse energy, Ej > 30 GeV and \rj\ < 2.5, excluding the transition region be- 
tween the barrel and endcap calorimeters, 1.44 < \ric\ < 1.57, where rjc is the pseudorapidity 
of the electromagnetic cluster. 

Signal events are required to have only one isolated lepton, whose origin is consistent with the 
reconstructed pp interaction vertex (Ts'l. Since the muon (electron) from a leptonic W decay 
is expected to be isolated from other high-px particles in the event, its track is required to be 
isolated from other activity in the event. This is done by requiring a relative isolation (/j-ei) 
less than 0.05 for muons and 0.10 for electrons. A looser cut is used in the electron case to 
allow for the increased amount of radiation close to the track. Relative isolation is defined as 
4ei = (^charged + ^neutral + ^photon) / Pt/ where pj is the transverse momentum of the lepton, and 
^charged/ ^neutral/ ^rid /photon are the sums of the transverse energies of the charged and neutral 
hadrons and the photons reconstructed in a cone of AR < 0.3 around the lepton direction, 
where AR = y/(A(/))2 + (A?/)2. The energy deposited by the lepton is explicitly removed from 
the sums by defining an exclusion cone of AR < 0.15 around the lepton direction. From studies 
based on Z decays in the data, the combined identification and isolation efficiencies for these 
selections are 83 ± 1% for muons and 75 ± 1% for electrons. 

Semileptonic tt events have at least four hadronic jets (from the hadronization of the bottom 
and light quarks). Charged and neutral hadrons, photons, and leptons are reconstructed using 
the CMS particle-flow algorithm [16J before they are clustered to form jets using the anti-A;T jet 
algorithm [17J with a cone size AR = 0.5. The jet clustering software used is FASTJET version 
2.4.2 |[l8l|l9l. At least one jet candidate with pj > 25 GeV and |//| < 2.4 is required, which 
must not overlap with a muon or electron candidate within AR < 0.3. Relative and absolute jet 
energy corrections |20| are applied to the raw jet momenta to establish a uniform jet response in 
Pj and ?] (with uncertainties of 3-6%, dependent on pj and rj). An offset correction, determined 
from simulation, is made to correct for the effects of event pileup (with an uncertainty < 1.4%). 
There is an additional b-jet energy scale uncertainty which accounts for the differences in re- 
sponse between Pythia and Herwig simulations (~3%). Lastly, a further 1.5% uncertainty 
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is added to allow for residual calibration differences during different run periods for a total of 
4.7-7.0%. 

The neutrino from the leptonic W decay escapes detection. Its presence is inferred from a 
sizeable transverse energy imbalance in the detector. The missing transverse energy (£t) is 
defined as the negative of the vector sum of the transverse energies (Ej) of all of the particles 
found by the particle-flow algorithm. This is used as an event selection variable in both the 
muon and electron analyses to suppress the background from QCD multijet events, and is 
required to be greater than 20 GeV in both channels. 

Because of the combined effects of the long b-quark lifetime (~1.5 ps) and the fact that the b 
quarks are produced with a significant boost, the decays of b-flavored hadrons are quite dif- 
ferent from the shorter lived hadronic states. Rather than having an origin consistent with the 
primary collision vertex, they can travel a measurable distance before decaying, resulting in a 
displaced decay vertex. The origin of the particles from the b decay is thus typically inconsis- 
tent with the primary vertex position. In the case of a semileptonic b-hadron decay, this results 
in the production of a lepton with a displaced origin, that is also embedded inside a jet. These 
characteristics can be used to identify b-quark jets and distinguish them from their non-b coun- 
terparts. Here we use a displaced secondary vertex algorithm to tag b decays and suppress the 
background from W(Z)+jet and QCD multijet events. The algorithm is described in detail in 
||2T1 and has a b-tag efficiency of 55% with a light parton (u, d, s, g) mistag rate of 1.5% for jets 
with Pj > 30 GeV in simulated QCD events. The event selection requires that at least one of 
the selected jets is b tagged. 

3 Signal and Background Modeling 

The efficiency for selecting lepton+jets signal events and the corresponding kinematic distribu- 
tions are modeled using a simulated tt event sample. The simulation is performed using Mad- 
Graph [22J, where the events containing top-quark pairs are generated accompanied by up 
to three extra partons in the matrix-element calculation. The parton configurations generated 
by MadGraph are processed with Pythia 6.4 l|23| to provide fragmentation of the generated 
partons. The shower matching is done using the Kt-MLM prescription [221. The generated 
events are then passed through the full CMS detector simulation based on Geant4 [24]. 

The production of W(Z) + jets events, where the vector boson decays leptonically has a similar 
signature and constitutes the main background. These are also simulated using MadGraph. 

We use a dynamical mass scale (Q^ scale) of (wJw/z)^ + (EPj')^ renormalization and 

factorization scales for both the W+jets and Z+jets simulations, and these scales are varied by 
factors of 2.0 and 0.5 in systematic studies. 

In addition to the Monte Carlo (MC) generation using MadGraph, QCD multijet samples 
were produced using Pythia. 

The QCD predictions for the top-quark pair production cross section are 157^24 pb in NLO [25] 
and 164+{° pb f9,^ or pb [7J in approximate NNLO. In each case the quoted uncertain- 

ties include the renormalization and factorization scale uncertainties, the uncertainties from 
the choice of parton distribution functions (PDF)s, and the uncertainty of the strong coupling 
constant (ag). For the scale uncertainty we have varied the renormalization and factorization 

scales by factors of 2 and 0.5 around the central choice of (2mt)^ + (X^ Pj')^ with nit = 172.5 GeV. 
The PDF and as uncertainties were determined by following the results obtained by using the 
MSTW2008 1261, CTEQ6.6 EZI, and NNPDF2.0 HSH sets and combining the results using the 
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PDF4LHC prescriptions ||29|. 

Single top quark production is described in terms of the cross sections in the s, t and tW chan- 
nels. The largest cross section is in the t channel where the predicted NLO cross section is 
at = 64.6+^;2pb from mcfm ||25l I30ti32l. The result is given for a scale of (mf)^ + {ZPt^)^ 
and an uncertainty that is defined in the same way as for top-quark pair production. Simi- 
larly, the cross sections for the tW and s channels are predicted to be cr^^ = 10.6 ± 0.8 pb and 
(Tg = 4,2 ± 0.2 pb 1 30 1, respectively. A measurement of the i-channel cross section has been per- 
formed by CMS and is found to be consistent with the value predicted by the standard model 
11331. 

The inclusive NNLO cross section of the production of W bosons decaying into leptons has 
been determined as c^w^/i/ = 31.3 ± 1.6 nb using FEWZ [34] (corresponding to a A:-factor of 1.3), 
setting renormalization and factorization scales to with mw = 80.398 GeV. The uncertainty 
was determined in a similar way as for top-quark pair production. Finally, the Drell-Yan pro- 
duction cross section at NNLO has been calculated using FEWZ as c^z/7*^/;("^;/ > 50 GeV) = 

3.0 ± 0.1 nb where the scales were set to m| with niz = 91.1876 GeV. 

The cross sections discussed above are used to normalize the simulated event samples. The 
profile likelihood fit yields corrections to these values. The normalizations for the final com- 
parisons with the data are determined from fits to the data in control samples. 

4 Cross Section Measurements 

The following subsections discuss the measurement procedure and the results obtained from 
the analysis of the muon, electron and combined channels. 

4.1 Fit Procedure 

To extract the tf cross section we perform a maximum likelihood fit to the number of recon- 
structed jets = 1-4, >5), the number of b-tagged jets (z = 1, > 2), and the secondary vertex 
mass distribution in the data. The secondary vertex mass is defined as the mass of the sum of 
the four-vectors of the tracks associated to the secondary vertex, assuming that each particle 
has the pion mass. It gives a good discrimination between the contributions from light and 
heavy flavor quark production (Fig. [l]). We fit the data to the sum of signal and background 
shapes using a binned Poisson likelihood. 

The templates for the fit are normalized to the expected event yields for 36 pb^^. The expected 
yield for each component in a given jet-multiplicity bin j and tag-multiplicity bin i is a function 
of the cross section, the jet energy scale (JES), the b-tag efficiency, and the background normal- 
ization. The JES and the Q^-scales affect the jet multiplicity, and the b-tag efficiency impacts 
the number of b tags. The background model depends on each of these and has additional con- 
tributions from the normalization and extrapolation from the sidebands in the data. Because 
these effects are expected to cause the largest uncertainties on the tt cross section measurement, 
and they are correlated, they are treated as nuisance parameters in the profile likelihood fit. The 
minimization of the likelihood thus provides simultaneous measurements of each contribution 
and the tt cross section. 

We determine the background normalizations from data sideband regions and extrapolate to 
the signal region using various models. In particular, the four major backgrounds are W+jets, 
Z+jets, single top quark, and QCD multijet production. 
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Figure 1: Secondary vertex mass distribution for bottom, charm, and light flavor jets. The 
bottom and charm templates are taken from simulated tt events and the light flavor shape is 
taken from simulated W+jet events. 

The templates for the W and Z contributions are normalized such that the NNLO predictions 
are equal to unity. During the profile likelihood maximization, the normalizations of these 
components are extracted. A full detector unfolding is not done, so this is not a meaningful 
measurement of the W(Z)+jets cross section, however the impact of the renormalization and 
factorization scales on the tt cross section is found to be smaller than that predicted by an ad 
hoc variation of the scales. 

The W and Z+ jets backgrounds come from V+h jets, V+c jets, and V+light flavor events. 
The same fc-factor of 1.30 (defined in Section |3]l is applied to all three flavors as inputs to the 
likelihood fit (although the three components are allowed to float independently in the fit). An 
additional electroweak background is single-top (s and t and tW channels) events. 

The shape of the jet multiplicity distribution (Njet) depends on the choice of the jet pj thresh- 
old, and thus is also sensitive to the jet energy scale (JES). In this sense, the fit is intrinsically 
able to determine the JES from the variations of the Njet distribution as a function of JES. The 
uncertainty in the b-tag efficiency is also extracted directly from the fit, by using the changes in 
the relative rates of 1-tag, and 2-tag events. A larger b-tag efficiency will result in events mov- 
ing from 1-tag to 2-tag samples. In contrast, an overall increase in all tag bins together would 
indicate an increase in the tt cross section. The combined in-situ measurement of the yields of 
principal backgrounds and parameters describing main systematic uncertainties leads to a sig- 
nificant improvement over analyses which use more conventional techniques for tagged cross 
section measurements. 

There are several "nonprompt-W" or "QCD" backgrounds for the muon and electron analyses. 
The QCD background in the muon-plus-jets channel comes from multijet events with heavy 
flavor decays, kaon and pion decays in flight, and hadronic punch-through in the muon system. 
Because these are difficult to calculate to the required precision we derive these backgrounds 
from the data. The normalization is determined by using a comparison of data and simulation 
in the data sideband region with < 20 GeV. The ratio is used to scale the predicted yields 
for £t > 20 GeV. The shapes and normalizations of the tt, W and Z+jets are well described 
by the Monte Carlo simulation and are modeled that way. The shape of the QCD component 
is derived from the nonisolated (Jrei > 0.2) data. Because of the correlations between and 
isolation, the templates for the QCD estimate from the data that are taken from nonisolated 
samples are modified using the shape taken from the QCD simulation. This treatment is similar 
in spirit to the QCD treatment in the recent CMS W and Z cross-section measurements [,35 J. 
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The QCD rate in each jet bin is constrained to the average of the true £t distribution of the 
nonisolated region, and the "modified" £t distribution after accounting for correlations from 
the Monte Carlo simulation. The QCD component in the fit is constrained to 100 percent of the 
rate or half the difference between the results, whichever is greater. 

The secondary vertex mass shapes for the fit are taken from the nonisolated data. Because of 
limited statistics, the > 3 jet sample is included as a single template, with separate normaliza- 
tions for each jet bin. 

Similarly, the electron-plus-jets channel is contaminated by photon conversions, jets with a high 
electromagnetic fraction, and heavy flavor decays. The normalization of this background is also 
estimated from a fit to the £t spectrum. However, the nonisolated sidebands do not accurately 
represent the shape of the distribution, and so the shape is determined by reversing at least 
two out of six of the electron identification criteria. The ensemble of these "marginal failures" 
makes a good representation of the shape of the £t background for this background source. 
This shape is then fit with the same £t procedure as in the muon case. 

There are alternative control samples for both the jet energy scale, and the b-tag efficiency. The 
jet energy scale is measured as described in Section |2] The uncertainty measured there is used 
in a Gaussian constraint on the jet energy scale in the likelihood (approximately 4%). The b- 
tagging efficiency and light quark mistag rate have been measured in an independent sample 
of QCD dijet events [36 1, and these values are input as Gaussian constraints on the parameters 
in the likelihood. Specifically, the b-tag efficiency scale factor is constrained to 1.0 ± 0.2, and 
the mistag rate scale factor is constrained to 1.0 ± 0.1. The technical implementation of these 
efficiencies in the likelihood is to weight the tagged jets in the simulation up or down by the 
data-to-simulation scale factor, and weight untagged jets in the simulation with zero weight. 

The number of predicted events with respect to each contribution is given by Eqns. for the 
tt signal and two of the W+jets backgrounds (W+b-jets and W+ light flavor). There are similar 
terms for the other W+jets backgrounds, the single top, and the QCD multijet production. Thus 
we have 



''Hhi,Rk tag) ■ P"^"'^S(wVRmistag) " P^"^ ( RjES ) 

N£(M-)=i<^wbb-Nmb(W> (2) 
P'' ''Hhi.Ri, tag) • P'"''''Hhi.Rr.ist.g) ■ PJ^S(/- ;,RjEs) ■ P^\hhRQ^) 

pmistag K^i3t,g) . pJES ^^^g) . pQ^ ) 



where K^i is the fitted scale factor for the NLO prediction for tt; i and / run over tags and jets, 
respectively; Ky^^ is the fitted scale factor for the NNLO prediction for Wbb (etc); N^^{i,j) 
is the number of events expected for sample X, derived from Monte Carlo and corrected with 
data-to-Monte-Carlo scale factors. The (/,;', Rx) factors are multiplicative functions account- 
ing for the relative differences with respect to the input expected yield, as a function of the 
assumed value Rx of nuisance parameter X (i.e., b-tag efficiency, jet energy scale, etc). These 
are interpolated from various configurations in the simulation with polynomials. The conven- 
tion chosen is that the nominal event yield is at Rx = (i.e., no variation in parameter X), and 
P{i,i,Rx) = 1-0 (i.e., multiplicative factor of 1.0 by default). The "+lcr" variation is at Rx = 1, 
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and the "—la" variation is at Rx = — 1- 

The fit minimizes the negative log likelihood, summing over the histogram bins (k) of the sec- 
ondary vertex mass, the number of jets (j), and the number of tags (z). The various constraints 
(described above) are included as Gaussian penalty terms on the variables, which are repre- 
sented by Cx- The full profile likelihood expression is 

{tagjetbins . . -i constraints _ ^ \2 1 

L L (ln^(Nf ^(z,;),Nr(M-))) - I L ^ \2 (4) 
i,i k I Cx ) 

where V is a Poisson probability that the predicted yield (N^"?) given by the various compo- 
nents statistically overlaps with the data (N*^) in each tag/jet bin given by 

lnP(;c,i/) = xlni/-i/-lnr(x + l) (5) 
where T{x) is the Gamma function. 

Table [ijshows a summary of all of the inputs to the profile likelihood, as well as the constraints. 



Table 1: Inputs to the profile likelihood, along with constraints. 



Quantity 


Constraint (%) 


b-tag Efficiency Scale Factor 

b-tag Mistag Scale Factor 

Jet energy scale relative to nominal 

W+jets renormalization/ factorization scales 

W+jets background normalization 

QCD background normalization 

Single-top background normalization 

Z+jets background normalization 


20 
10 
4 

+100 
-50 

unconstrained 
100 
30 
30 



There are also a number of systematic uncertainties that are not included directly in the profile 
likelihood and hence are taken as additional systematic uncertainties outside of the fit result. 
The largest of these is the systematic uncertainty due to the overall luminosity determination. It 
has also been shown on independent samples of Z — > ee and jiji events, that the efficiencies 
for triggering, reconstructing, and identifying isolated leptons of this type are very similar in 
the data and simulations. We have corrected for the small differences observed. The effect 
of these uncertainties are not included in the profile likelihood, and hence are taken as an 
additional systematic uncertainty of 3%. 

There are a number of theoretical uncertainties in the signal modeling that are not included in 
the profile likelihood at this time. They include differences in the tt signal due to renormaliza- 
tion and factorization scales, the amount of initial and final state radiation present, the parton 
distribution function model, and the matching scale for the matrix-element to parton-shower 
matching scheme. These are computed from dedicated simulated samples by varying the the- 
oretical parameters of interest according to conservative variations around the reference value. 
The exception is the parton distribution functions, which are varied by reweighting the sample 
according to variations in the underlying parton distribution function parameterizations [29]. 
The numerical impact of each of these is taken as a systematic uncertainty. Specifically these 
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are 2% for tt modeling; 2% for initial and final state radiation modeling; 1% for the matrix- 
element to parton-shower matching in the tt; and 3% for the parton distribution function dif- 
ferences. 

In all of the cases that are described below, the robustness of the statistical procedure is demon- 
strated with a priori pseudo-experiments where the expected yields and the parameters in the 
profile likelihood are sampled randomly according to Poisson or Gaussian statistics (as ap- 
propriate). In cases where the true frequentist statistical coverage is not achieved due to the 
limitations of the profile likelihood method, coverage is assured by correcting for the slight 
biases (of order 1-2% in the central values and/ or uncertainties). 

4.2 Muon + Jet Analysis 

For the muon channel we make two modifications to the basic event selection discussed above. 
Instead of using the isolation cut of Jj-ei < 0.05, we require that the selected muon tracks are 
found to be isolated by the particle-flow reconstruction and pass a cut of Jrei < 0.15. The results 
of the fit are shown in Table |2j The fit considers events with one tag (1-tag) and two and 
more tags (2-tag) separately, giving nine jet-tag loins' (subsamples) which are fit by the joint 
likelihood. Table|2]lists the observed and fitted rates for each jet-tag bin. 



Table 2: Results of the fit for muon+jets events with at least 1 b tag. 





Data 


Fit 


tt 


t(t) 


W 


W 


W 


Z 


QCD 












+b jets 


+c jets 


+qq 


+jets 




1 jet 1 b tag 


505 


504.0 


13.3 


25.0 


94.2 


255.1 


81.9 


13.9 


20.6 


2 jets 1 b tag 


314 


318.2 


51.0 


29.4 


82.6 


97.7 


35.0 


7.3 


15.1 


3 jets 1 b tag 


166 


158.5 


78.3 


14.8 


29.5 


21.9 


10.4 


2.8 


0.8 


4 jets 1 b tag 


85 


89.2 


60.6 


4.9 


12.8 


5.5 


3.3 


1.3 


0.8 


> 5 jets 1 b tag 


45 


43.8 


34.6 


1.5 


4.6 


1.8 


0.9 


0.2 


0.2 


2 jets >2 b tags 


29 


24.1 


14.7 


3.3 


5.3 


0.5 


0.2 


0.2 


0.0 


3 jets >2 b tags 


37 


44.0 


35.2 


3.8 


3.9 


0.8 


0.0 


0.2 


0.0 


4 jets >2 b tags 


41 


41.0 


36.2 


1.9 


2.4 


0.4 


0.0 


0.1 


0.0 


>5 jets >2 b tags 


27 


26.0 


24.0 


0.8 


1.0 


0.1 


0.1 


0.1 


0.0 


Total 


1249 


1248.8 


347.8 


85.4 


236.2 


383.8 


131.8 


26.1 


37.6 



Table |3] lists the systematic uncertainties from the fit. These include both the theoretical uncer- 
tainties from the tt modeling and the corrections used to match the simulations to the data and 
give a total uncertainty of 4.3% for the tt signal model. The unclustered energy in the detector 
results in an additional resolution uncertainty of <1% on the scale. We combine these with 
the data-simulation uncertainties due to the jet energy scale and jet resolution modeling, the 
b-tag efficiency and mistag rate and obtain a total systematic uncertainty of 12.5%. For illustra- 
tive purposes, in Table |3} we have broken up the pieces of the profile likelihood and quote the 
uncertainties due to the individual contributions. These are the result of fixing all of the other 
parameters of the likelihood and only allowing the chosen term to vary. 

This yields a cross section measurement of 

cT-tf = 145 ± 12 (stat.) ± 18 (syst.) ± 6 (lumi.) pb, (6) 

where the last uncertainty corresponds to the 4% uncertainty on the total integrated luminosity 
(BJ. The fit provides in-situ measurements of the scale factors for both b tagging and the jet 
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Table 3: List of systematic uncertainties for the muon + jet, electron + jet, and combined 
analyses. Due to the correlation between parameters in the fit, the combined number is not the 
sum of the squares of the contributions. 



Source 


Muon 


Electron 


Combined 




Analysis 


Analysis 


Analysis 


Quantity 


Uncertainty 


(%) 


Lepton ID/ reco/ trigger 




3 




£t resolution due to unclustered energy 




< 1 




tt+jets scale 




2 




ISR/FSR 




2 




ME to PS matching 




2 




PDF 




3 




Profile Likelihood Parameter 


Uncertainty (%) 


Jet energy scale and resolution 


10 


9 


7 
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9 


8 


8 


W+jets scale 


4 


3 


9 


Combined 


13 


12 


12 



energy scale. We obtain a result of 98 ± 6% for the b-tag scale factor and 92 ± 9% for the jet 
energy calibration. The scale factors for the W+b jets and W+c jets components indicate that the 
contributions in the data may be larger than what is expected by the predictions. For the W+b 
jets component we find a cross section scale factor of 2.6 and for the W+c jets contribution 
we obtain 1.3 ^q 2. It is also found that the W+jets data are slightly harder than the central value 
of the renormalization and factorization scales chosen. 



4.3 Electron + Jets Analysis 

The analysis in the electron channel is performed in the same way as for the muon case. The 
results are shown in Table |4] The fit was performed in the same manner as for the muon 
channel, resulting in nine jet-tag 'bins' (subsamples), which were fit by a joint likelihood, as 
described in Section 4.1 Table |4] lists the observed and fitted rates for each jet-tag bin. Note 
that the fit parameters are unbounded to avoid problems with bias. This can result in negative 
values for the event counts, as in the case of the 4 jets 1 tag yield for QCD in Table |4] 



Table 4: Results of the electron+jets fit for events with at least 1 b tag. 





Data 


Fit 


tt 


t(t) 


W 
+b jets 


W 
+c jets 


W 
+qq 


Z 
+jets 


QCD 


1 jet 1 b tag 


388 


389.8 


6.0 


14.1 


42.4 


249.5 


26.9 


3.1 


47.9 


2 jets 1 b tag 


252 


245.3 


31.7 


21.0 


44.4 


104.0 


14.8 


3.3 


26.2 


3 jets 1 b tag 


159 


156.0 


62.7 


12.2 


23.5 


34.4 


5.1 


2.1 


16.0 


4 jets 1 b tag 


71 


80.7 


60.6 


4.8 


8.2 


9.2 


1.4 


0.8 


-4.3 


>5 jets 1 b tag 


57 


52.1 


40.9 


1.6 


2.6 


3.0 


0.5 


0.4 


3.0 


2 jets >2b tags 


14 


19.9 


9.4 


2.3 


4.7 


1.2 


0.1 


0.2 


2.0 


3 jets >2 b tags 


39 


38.1 


29.1 


3.1 


3.8 


0.9 


0.0 


0.2 


0.9 


4 jets >2b tags 


37 


41.3 


37.1 


1.9 


1.9 


0.4 


0.0 


0.1 


0.0 


>5 jets >2 b tags 


37 


30.7 


28.8 


0.9 


0.8 


0.2 


0.0 


0.1 


0.0 


Total 


1054 


1053.8 


306.3 


61.8 


132.1 


402.7 


48.9 


10.3 


91.6 
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4 Cross Section Measurements 



The resulting cross section is 



o-ft = 158 ± 14 (stat.) ± 19 (syst.) ± 6 (lumi.) pb. 



(7) 



From the fit we obtain a result of 97 ± 6% for the b-tag scale factor and 103 ± 8% for the jet 
energy calibration. The scale factors for the W+b jets and W+c jets components are 1.4 ^ 
and 1.4 t^^'^, respectively. These are in agreement with the results from the muon channel. The 
contributions to the systematic uncertainty are summarized in Table |3] 

4.4 Simultaneous Muon and Electron Channel Analysis 

Having established the consistency of the separate channel measurements, we now proceed to 
perform a combined fit to both channels. To establish our best measurement, we repeat the 
fit procedure and apply it simultaneously to the data in both the electron and muon channels. 
We find that the resulting fitted event yields in each tag category are in good agreement with 
those obtained from the separate channel fits (Tables |2] and |4]). Figure |2] shows the comparison 
of the corresponding observed and fitted vertex mass distributions. Figure |3] shows the data 
for > 3 jets and > 1 b-tag, and the fit results for the total transverse energy of the event (Hj), 
the missing transverse energy (E^^^), and the transverse mass of the W (M^). We find good 
agreement in all cases. 

. Data Jttjsingle Top |w^-jets Jw4c-jets [^W4iF-jets Jz-ljets [^OCD 



CMS, VS =7TeV 36 * pb of Muon Data 




^250= 
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Figure 2: Results of the combined muon and electron channel fit. The muon channel is shown 
on the left and the electron channel on the right. The plots on the top are for exactly 1 b tag 
and those on the bottom are for >2 b tags. The histograms within the top panel correspond 
to events with 1, 2, 3, 4 and >5 jets, respectively, while the bottom panel shows histograms 
corresponding to events with 2, 3, 4 and >5jets. 

The correlation matrix for the combined fit is listed in Table H] All of the terms are as defined 
in the text. The combined analysis cross section measurement is 
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Table 5: Correlation matrix of the fit to the combined electron and muon data samples. 









w 

vv 


w 

vv 


w 

vv 


7 


\d 


b tag 


J i-jj 


R , 

^^mistag 








+b jets 


+c jets 


+qq 


+jets 










tt 


1.0 


-0.1 


0.1 


0.5 


0.2 


0.0 


0.3 


-0.7 


-0.6 


0.0 


t(t) 


-0.1 


1.0 


-0.3 


0.0 


0.0 


0.0 


0.1 


0.1 


0.1 


0.0 


W+b jets 


0.1 


-0.3 


1.0 


0.4 


0.4 


0.0 


0.6 


-0.2 


0.0 


0.0 


W+c jets 


0.5 


0.0 


0.4 


1.0 


-0.1 


0.0 


0.6 


-0.5 


-0.2 


0.0 


W+qq 


0.2 


0.0 


0.4 


-0.1 


1.0 


-0.1 


0.5 


-0.2 


-0.2 


-0.3 


Z+jets 


0.0 


0.0 


0.0 


0.0 


-0.1 


1.0 


0.0 


0.0 


0.0 


0.0 


Q' 


0.3 


0.1 


0.6 


0.6 


0.5 


0.0 


1.0 


-0.2 


0.1 


0.0 


btag 


-0.7 


-0.1 


-0.2 


-0.5 


-0.2 


0.0 


-0.2 


1.0 


0.3 


0.0 


JES 


-0.6 


0.1 


0.0 


-0.2 


-0.2 


0.0 


0.1 


0.3 


1.0 


0.0 


^mistag 


0.0 


0.0 


0.0 


0.0 


-0.3 


0.0 


0.0 


0.0 


0.0 


1.0 



(7tt = 150 ± 9 (stat.) ± 17 (syst.) ± 6 (lumi.) pb, (8) 

which is in good agreement with both the separate channel measurements and those from the 
cross-check analyses discussed below. The corresponding summary of the systematic uncer- 
tainties is given in Table |3j We obtain a result of 96 for the b-tag scale factor which agrees 
well with the result obtamed in [21]. For the jet energy scale we obtain a result of 107 ± 6% 
indicating that the data may prefer a small increase in the jet energy calibration. The scale fac- 
tors for the W+b jets and W+c jets components indicate that the contributions to the data may 
be larger than what is predicted by the scaled NNLO predictions of 110 pb and 3.0 nb. For the 
W+b jets contribution we find a cross section scale factor of 1.9 5, which is similar to recent 
observations at the Tevatron ||37f|39l . The result for the W+c jets contribution is 1.4 ± 0.2. 

5 Cross-check Analyses 

As a cross-check of our results, we have performed a series of independent analyses in both 
the muon and electron channels. A summary of these is given in the following subsections. 
These use not only different analysis techniques but also different methods to suppress the 
backgrounds from W + jets and QCD multijet events. Each analysis requires at least three jets 
with Pj > 25 GeV and has no requirement on the amount of missing transverse energy. In 
addition, two different tagging algorithms are used for the two analyses in the muon chan- 
nel. The Neural Network analysis uses a track counting algorithm which counts the number of 
tracks nonassociable to the primary vertex l|2T| and the second analysis uses the muons from 
semi-leptonic b decays to tag b jets. Each analysis has significantly different systematic uncer- 
tainties from the analysis presented above and thus provides a good test of the robustness of 
the measurements. While there are some differences in the selected event samples, we do not 
attempt to combine the results because of the substantial overlap. 

5.1 Neural Network Analysis with a Track-counting b-tagger 

The first cross-check is performed in the muon channel. It makes use of a multi-layer per- 
ceptron neural network to distinguish tt signal events from the backgrounds after requiring a 
muon and three jets to pass the selection criteria. The network discriminant is built from the 
analysis of three input variables: the pseudorapidity | j/^ | of the muon, the distance AR12 in rj-cp 
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Figure 3: Kinematic distributions of the total transverse energy (Hj), the missing transverse 
energy (E^^®) and the transverse mass of the W (£t). The plots are for > 3 jets and > lb tag 
and the comparison histograms correspond to the fitted values. The muon charmel is shown 
on the left and the electron channel is on the right. 

space between the two highest-px jets in the event, and a boolean variable indicating the pres- 
ence of at least one b-tagged jet. The network was trained using simulated samples of signal 
and background events and the cross section was determined by fitting the sum of signal and 
background templates to the data. 

The signal discriminant was generated from a simulated tt sample which was corrected to 
match the jet energy resolution observed in the data. An additional flavor-dependent correction 
was applied to the simulated jets to account for differences in b-tag efficiencies between data 
and simulation. The template shapes for QCD, W+jets, and Z+jets were produced directly 
from control regions in the data. These were chosen using variables that are only very loosely 
correlated to the network parameters so that they do not bias the operation of the network. 

Because the event topologies of tt and single-top events are so similar, there is very little differ- 
ence in the shapes of the discriminant for the two samples. In order to avoid pathological fit 



5.2 Muon Channel Analysis using a Muon-in-Jet b-tagger 
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CMS,\/s = 7 TeV, 36 pb"^ of Muon Data 
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Figure 4: Results of the template fit to the neural network discriminant for 36 pb^ of data. 

yields due to this similarity, the number of single-top events in the fit was constrained to the 
expected yield from the NLO single-top cross section with a 30% uncertainty. The QCD and 
W/Z+jets likewise share similar shapes in their discriminants. This ambiguity was resolved by 
constraining the QCD fit yield according to its leading-order (LO) cross section, with a 100% 
uncertainty due to expected differences between the actual QCD cross section and its LO cal- 
culation. However, as the QCD contribution is small, it has little effect on the final result. 

The fit of the neural-network discriminant to events passing selection cuts in the data is shown 
in Fig. |4] A study using pseudo-experiments with simulated data indicates that the fitter in- 
troduces a —3.1% bias in the tt yield. This arises due to the use of control regions for the 
QCD, W+jets, and Z+jets templates. After correcting for this, we obtain a tt event yield of 
369 ± 36 (stat.) events. 

The systematic uncertainties are evaluated by shifting the simulation by each systematic un- 
certainty and re-evaluating the result. The resulting pseudo-data are fit with the nominal dis- 
criminator templates for tt and single top, and systematically shifted templates for QCD and 
W/Z+jets. The dominant uncertainty of 1|5% comes from the b-tag efficiencies. When com- 
bined with the other contributions, we find a combined systematic uncertainty of ]'^^8°^°- After 
combining this with the 4% uncertainty in the recorded integrated luminosity IHI, we obtain a 
tt cross section of 151 ±15 (stat.) (syst.) ±6 (lumi.) pb. 

5.2 Muon Channel Analysis using a Muon-in-Jet b-tagger 

We have also performed an analysis in which b jets are identified by the presence of a noniso- 
lated muon. The event selection requires at least three well reconstructed jets of which at least 
one contains a muon with pj > A GeV and \ f]\ < 2.4. The backgrounds from events which come 
from the decay of a J/t/^, Y, or Z boson, are excluded using selections on the muon pair invariant 
mass. 

The background from W+b jets and W+c jets with a semimuonic decay and W+light flavor 
events with a jet that is misidentified as containing a semimuonic decay are estimated from 
data. We calculate the track taggability from a sample of 7+jet events. This taggability is 



14 



5 Cross-check Analyses 



convoluted with the track distribution in the jets of the pre-tagged event sample to predict the 
number of tagged events that should arise from W+jets, Ntag,pred- We correct for the fraction 
of events in this sample that are due to QCD multijets, Fqcd, to obtain Nw+jets = Ntag,pred ■ 
(1 — Fqcd)- We estimate Fqcd using the method discussed below. We assign a 30% systematic 
uncertainty, based on studies of the taggability parameterization as applied to independent 
control samples, to the tag-rate prediction to account for the simulation uncertainties. 

The QCD background was estimated by calculating the fraction of QCD events in the pre- 
tag sample. Because of the enrichment in bb and cc events after requiring the primary muon, a 
correction factor k is applied to correct for this. The QCD background is then Nqcd = Ntag,pred " 
Fqcd ■ k. After taking into account the uncertainties in the calculation, we assign a systematic 
uncertainty in the QCD background of 60%. 

The background from Drell-Yan events that survive the Z veto is estimated from simulation. 
An estimate from data similar to the one described in reference f40] was employed to assign a 
systematic uncertainty to the prediction from simulation. This gives a systematic uncertainty 
of 17% on the Z mass veto correction. The remaining backgrounds from diboson and single- 
top production were estimated from simulated samples that were normalized to the theoretical 
NLO cross sections. Each of these is assigned a 30% systematic uncertainty. 

The selection efficiency for signal events prior to the muon-tag requirement is evaluated using 
simulated tt events and the method described previously. The combined systematic uncer- 
tainties on the acceptance, due to the simulation uncertainties, is 9%. Because of the rela- 
tively soft Pj spectrum of the tag muons, the reconstruction and identification efficiencies were 
checked using a tag-and-probe analysis f35l of ]/tp — ?► events. The results agree with 

the simulation to within 1%, so we take the efficiency for finding a tag muon and the mistag 
efficiency directly from simulation. This gives a tagging efficiency for tt events of 25.4 ± 0.1% 
excluding the resonance veto requirements, and 23.9 ± 0.1% including them. The systematic 
uncertainty on this efficiency is conservatively taken to be 10%. Before calculating the tt cross 
section, the number of predicted tagged events is corrected for the presence of tt events in the 
pre-tag sample. This correction is performed iteratively. Combining these effects gives a total 
systematic uncertainty on the tt cross section due to the background calculations of 17.9%. 

Figure |5] shows the observed jet multiplicity distribution, together with the signal and back- 
ground predictions. Here, the tt component is normalized to the NLO cross section II25II and 
the backgrounds are those computed above. The combined signal and background predictions 
are in good agreement with the data for each jet multiplicity. 

Combining the data yield, the background estimation, the tt acceptance, and the tagging ef- 
ficiency with all of the associated uncertainties, gives a cross section result of 163 ± 21 (stat.) 
± 35 (syst.) ± 7 (lumi.) pb. This measurement provides a valuable cross-check of the muon+jets 
result, as its systematic uncertainties are almost independent of those in the reference analysis. 

5.3 Electron Channel Cross-Check Analysis 

Our third cross-check analysis is a simple counting analysis using the same displaced vertex 
tagger as our reference analysis. The event selection requires the presence of at least three well 
reconstructed jets, of which one or more is required to be b tagged. We estimate the QCD 
background from the data using control regions in isolation and the W+jets background from 
a combination of data and simulation. To calculate the tf event yield, these are subtracted from 
the data, along with smaller contributions from single-top, Drell-Yan, and diboson production, 
which are estimated from simulation. 



5.3 Electron Channel Cross-Check Analysis 
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CMS,\/s = 7 TeV, 36 pb"' of Muon Data 




Jet Multiplicity 

Figure 5: Comparison of data and simulation for the jet multiplicity for events with > 1 muon- 
in-jet tag. The W+jets and QCD contributions are normalized to the data-driven predictions. 
The hatched area shows total uncertainty on the background prediction. 



The selection efficiency is calculated for tt and for W+jets, and the efficiencies to tag individ- 
ual jets are measured in bins of jet-px/ separated into h, c, and u, d, s, g (light) partons. The 
uncertainties from the data are used in systematic studies as well as the statistical uncertainty 
from the jet efficiency calculations. The resulting predicted b-tag rates applied to simulation of 
signal and backgrounds are in good agreement with the data (see Fig.|6]). 

CMS, is=7 TeV, 36 pb"^ of Electron Data 



CO 
-(— ' 


350 


CD 




> 




LU 


300 




250 




200 




150 




100 




50 








- e + >3 jets 






[ + 


-•- Data 




■ Single Top 

□ W^lv - 

■ Z/y*^lT 

□ QCD - 

□ VV : 




















1 2 


3 



Secondary vertex tag multiplicity 

Figure 6: Number of observed and predicted secondary vertex tags in events with > 3 jets. 



The QCD background estimation is performed using a fit to the electron I^^i distribution in 
the data above the standard selection where the QCD multijet events dominate. A Gaussian 
function is used for the central value as it best fits the data, and the fit uncertainty is estimated 
by varying the Jj-ei range used in the fit and by using alternate fit functions. The systematic 
uncertainty from tt contamination is evaluated by subtracting the number of tt events predicted 
to be inside the fit region using simulation. 

The background due to W+jets is estimated using a technique motivated by Berends-Giele 
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scaling l|4T]| to measure the jet multiplicity distribution before b-tagging. In strict Berends- 
Giele scaling, the ratio C(n) of events with W+ > n jets to events with W+ > (n + 1) jets is 
expected to be independent of n. To account for an observed deviation from perfect scaling, in 
this estimate, C(n) is extracted as a linear function of n, with slope taken from the simulated 
event sample for 1 < n < 3. The W+b jets content in the simulation is scaled up by a factor of 



two, as determined in Section 4.4 



A W+jets jet multiplicity spectrum in data is prepared by applying an additional requirement 
on the transverse missing energy, > 20 GeV, to suppress the QCD background. The number 
of QCD events remaining in these data are then subtracted by refitting the Jj-ei distribution as 
above. The remaining non-W backgrounds and the tt content are subtracted from the data 
sample using Monte Carlo predictions. The normalization of C(n) is then fitted to this data 
sample. The scale factors obtained are C(l) = 4.91 ± 0.13 and C(2) = 5.35 ± 0.16, where the 
uncertainties are statistical. Using C(n) the estimate of the W+jets content in the data with > 3 
jets before b tagging (AfP'-eteg'data) jg 



J. Tpretag,data 
J. Tpretag,data _ W+jets,>ljets 
-'^W+jets,>3jets ~ C(1)C(2) ' 

The number of tagged W+jet events is then found by applying the selection efficiencies de- 
scribed above. 

Systematic uncertainties for the W+jets background estimate are derived by repeating the above 
process varying tt content by 30% as well as varying the normalization of W+b jets by 50% and 
W+c jets content by factors of 2 and 0.5. Additional fits are also done to Monte Carlo samples 
where the factorization scale is doubled or halved, as a systematic study. The backgrounds 
from single-top, Drell-Yan, and diboson production are taken from simulation with a 30% sys- 
tematic uncertainty. 

The background-subtracted and efficiency-corrected yield gives a cross section measurement 
of 169 ± 13(stat.) (syst.) ± 7 (lumi.)pb. The dominant systematic uncertainties are the 
uncertainties on the b-tag scale factors, and the uncertainty on the jet energy scale. This analysis 
uses a significantly different analysis strategy from the the reference electron+jets analysis but 
obtains a very similar cross section. Thus it provides a good cross-check of the method and the 
results. 

5.4 Additional Cross-Checl( 

As a further cross-check, the tt cross section was also measured in the muon channel using 
a simple counting analysis which used the same vertex tagger as our reference analysis. The 
analysis is similar to that of the electron channel analysis described in the previous subsection. 
The data were selected using a relaxed l^^i cut of < 0.1 and a missing transverse energy require- 
ment of ?T > 20 GeV. The QCD multijet background was measured from pre-tagged data using 
the matrix method [42] and the W + jet background was derived from background-and-signal- 



subtracted data using the Berends-Giele method described in the Section 5.3 After requiring at 
least three well reconstructed jets with at least one jet with a b tag and applying a correction of 
factor of 2 ± 1 for the W + heavy flavor content of the simulated data, we find good agreement 
between the data and simulations. The resulting cross section is in good agreement with our 
reference result. 
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6 Combined CMS Measurement 

In addition to the results from this analysis, CMS has also performed a measurement in the 
dilepton decay channel Il43ll , where we measured a cross section of 

dileptons: cTjj = 168 ± 18 (stat.) ± 14 (syst.) ± 7 (lumi.) pb. 

To produce a final CMS result, we combine this result with the lepton + jets measurement using 
a profile likelihood method. This procedure uses the vertex mass templates information from 
this analysis and adds the dilepton measurements as single bin templates. For the dilepton 
channels, six statistically independent inputs are used so the uncertainty correlations can be 
handled correctly. These are the e}i events with > 2 jets and no b tag requirement, the ef/, ee, 
and mi events with only 1 jet, and the ee and ^^ events with > 2 jets and at least 1 b tag. 
Because different b taggers were used in the two analyses, the b-tag uncertainties are treated 
as uncorrelated. The 3% PDF uncertainty from the lepton + jets measurement is also treated 
as uncorrelated because it does not relate to anything in the dilepton analysis. The remaining 
systematic uncertainties (JES, lepton efficiency, renormalization and factorization scales, ME to 
PS matching, and ISR/FSR) are all assumed to be correlated. The method was verified by using 
pseudo-experiments, from which it was determined that the likelihood calculation resulted in a 
10% underestimate of the errors. We correct for the underestimate in the final result and obtain 
a combined measurement of 

CMS combined: a^^ = 154 ± 17 (stat. + syst.) ± 6 (lumi.) pb. 

Figure [7| shows the comparison of the separate and combined CMS IITOl |43| measurements of 
the production cross section. The inner error bars on the data points correspond to the statis- 
tical uncertainty, while the outer (thinner) error bars correspond to the quadratic sum of the 
statistical and systematic uncertainties. The outermost brackets correspond to the total uncer- 
tainty, including a luminosity uncertainty of 11% (4%) for the 3 (36) pb^^ results, respectively, 
which is also added in quadrature. 

Also shown are NLO and approximate NNLO QCD calculations, for comparison. These were 
computed using the HATHOR program |5| using the calculations from ||6|. A common factor- 
ization and renormalization scale of Q = mt = 173 GeV was used for the calculations together 
with the MRTSW 2008 NNLO (NLO) parton distribution functions. The scale uncertainty was 
determined by independently varying the two scales by factors of 2 and 0.5 and taking the max- 
imum variation as the uncertainty. The PDF uncertainty corresponds to the 90% confidence 
level (C.L.) uncertainties for the parton distribution functions p9l| . This is added in quadrature 
to the scale uncertainty. 

In Fig.|8]we compare these and the pp results from the Tevatron ['37','38','44','45'| to the theoretical 
predictions as a function of ^/s. These are from the approximate NNLO QCD calculations 
referred to above and the width of the error band corresponds to the variation resulting from 
changing the scale up and down by a factor of 2. We find good agreement between the data 
and the theory in all cases and we note that the combined measurement is already more precise 
than the NLO QCD prediction. 

7 Summary of Results 

We have presented the results of a new analysis of the tt production cross section at ^/s = 7 TeV 
using data recorded by the CMS detector at the LHC during 2010 corresponding to an inte- 
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Figure 7: Comparison of the CMS measurements for the tt production cross sections and the 
QCD predictions for a/s = 7 TeV 

grated luminosity of 36 pb . Using muon and electron+jets and using b tagging to suppress 
the backgrounds, we measure cross sections of 

jW+jets: o-^^ = 145 ± 12 (stat.) ± 18 (syst.) ± 6 (lumi.) pb, 
e+jets: a^^ = 158 ± 14 (stat.) ± 19 (syst.) ± 6 (lumi.) pb, 

from the separate channels. The combination of these gives a cross section of 

Z+jets: CTjj = 150 ± 9 (stat.) ± 17 (syst.) ± 6 (lumi.) pb. 

When combined with the CMS dilepton measurement, we obtain an improved cross section 
measurement of 

CMS combined: (r^i = 154 ± 17 (stat. + syst.) ± 6 (lumi.) pb. 

The measurements are in good agreement with the QCD predictions of 164^^^ pb fS*. ^ and 
163+JJpb [7] which are based on the full NLO matrix elements and the resummation of the 
leading and next-to-leading soft logarithms. 
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